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We analyze the effects of a nonlinear cubic perturbation on the (5-Kicked Rotor. We consider 
two different models, in which the nonlinear term acts either in the position or in the momentum 
representation. We numerically investigate the modifications induced by the nonlinearity in the 
quantum transport in both localized and resonant regimes and a comparison between the results for 
the two models is presented. Analyzing the momentum distributions and the increase of the mean 
square momentum, we find that the quantum resonances asymptotically are very stable with respect 
to the nonlinear perturbation of the rotor's phase evolution. For an intermittent time regime, the 
nonlinearity even enhances the resonant quantum transport, leading to superballistic motion. 

PACS numbers: 05.45.-a,03.65.Ta,42.50.Vk 



I. INTRODUCTION 

Recent and ongoing experiments I'] have started to in- 
vestigate the interplay between the many-body induced 
self-interaction in an ultracold atomic gas and an exter- 
nal driving induced by time-dependent optical potentials. 
The natural setup is to use a Bose-Einstein condensate 
of alkali atoms, where the nonlinearity parameter can be 
tuned 0, 13 , and pulsed optical lattices can be used to 
impart momentum kicks on the atoms. For such a setup, 
the Gross-Pitaevskii (GP) equation 2,'3'| provides a good 
description of the system, as long as the nonlinearity is 
not too large, as a study of the stability of linearized 
excitations around the GP solution has shown 0, 0| . 

In this paper, we analyze the evolution of a cubic non- 
linear Schrodinger equation, as present in the GP model, 
under the perturbation of time-periodic (5-kicks 
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where and n — — i^ are the position and the con- 
jugated momentum of the system; we chose units such 
that h = 1 and the motion is considered on a ring with 
periodic boundary conditions ■)/)(?? -f27r) = The pa- 

rameters u and k are the nonlinearity coupling and the 
kicking strength, respectively. 

In atom optics experiments, the (5-kicked rotor has 
been realized with an ensemble of laser-cooled, cold 
atoms 0, or recently also with an ultracold Bose- 
Einstein condensate ^ IJ , periodically driven with a stand- 
ing wave of laser light. With the wave number of the laser 
/cl, the experimental variables are easily expressed in our 
units by noting that momentum is usually measured in 
two photon recoils {2hkL), and position in units of the in- 
verse wave number of the standing wave (l/2fci). Hence, 
the scaled variables n and the physical ones are 



related by = 2kL^' and n = p' /2kLh 0, S 0. Our 
choice of units makes all the relevant quantities (includ- 
ing the ones plotted in figures) dimensionless. 

Owing to the periodicity of (5-kick perturbation, the 
time t is measured in number of periods r and the evo- 
lution of the wave function of the system over a time 
interval r is described by the operator U (r) . 

A cubic modification of linear Schrodinger dynamics 
for the (5-kicked rotor may be accomplished by two dif- 
ferent models, both considered in the present paper. The 
correct way to approximate the evolution of the nonlinear 
Schrodinger equation is to evaluate the nonlinear term in 
the position representation In following we will 
refer to this first model as model 1 (Ml). 
Since the Hamiltonian operator presents the time depen- 
dent nonlinear part u\il'\'^ipj iu ths numerical integration 
of Eq. , the lowest order split method is used and 
U is approximated by the time ordered product of evolu- 
tion operators (Trotter-Kato discretization 12]) on small 
time steps t/L (with L integer): 
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In the numerical simulations, we use a finite Fourier ba- 
sis of dimension N: the discrete momentum eigenval- 
ues lie on the lattice p = {m — y) ^ud the continuous 
angle variable is approximated by = ^(m — 1) with 
TO G Z, 1 ^ rn ^ iV. Shifting between the coordinate 
and momentum representations, in the evaluation of the 
operator W^'^It), requires 2L Fast Fourier transforms of 
A^-dimensional vectors for each kick. In order to get sta- 
ble numerical results, the splitting interval t/L has to be 
reduced when increasing the nonlinear coupling constant 
u; typical values of the number of steps per period range 
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between L = 80000 and L = 5000000. Therefore, inves- 
tigating either the effects of strong nonlinearities or the 
dynamics of the system over long times is computation- 
ally quite expensive with (Ml). 

A second model (known as the Kicked Nonlinear Rota- 
tor and called in the following model 2 (M2)) was intro- 
duced in [l3j . This model, being a much simpler variant 
of the Kicked Rotator Model (KR) JJJ, allows to per- 
form faster and more efficient numerical computations. 
The evolution operator over one period r for (M2) is 
given by 

_ g-ifecosWg-.ri^g.«r|^„|= ^3) 

where ipn indicates the n-th component of the wave func- 
tion of the system in the momentum representation. The 
change in the phase of each component V'n of the state 
vector, introduced by the nonlinear term between two 
kicks, is proportional to the amplitude of the component. 
In (Ml) instead, the phase acquired at each instant by 
the wave function involves all the Fourier components 
and the phase factor has as the n-th Fourier component 
5^ V'm-i-nV'm- The two modcls coincide only if the 
wave function of the system is a plane wave of fixed mo- 
mentum; in this case, the relation between the nonlinear 
coupling constants in (Ml) and (M2) is u = u/2n. 

Both models (Ml) and (M2) are nonlinear generaliza- 
tions of the KR and reduce to the KR in the limit u ^ 0. 
Depending on the commensurability of the period t of 
the (5-kick perturbation with 47r, the KR displays differ- 
ent regimes, deeply studied both theoretically pi^llSlllq 
and experimentally 0,0- The quantum resonant regime, 
corresponding to values of r being rational multiples of 
47r, is characterized by a ballistic transport: the mean 
energy of the system grows according to a parabolic law 
[isl fl^. For generic irrational values of r, the average 
energy grows linearly in time only within a characteris- 
tic time (break time) , after which dynamical localization 
sets in and the diffusion is suppressed p^ . 

In this paper, we analyze in detail how the presence of 
the nonlinearity affects the general properties of trans- 
port in regimes that correspond to the localized (Section 
II, which is essentially a warm up exercise in which the 
results of [l^ Eire reproduced and new numerical results 
about prefactor scaling is presented) and resonant (Sec- 
tion III) one of the KR. We focus the attention on the 
growth exponent of the mean energy and on how the 
diffusion coefficient or the rate of ballistic transport de- 
pend on the strength of the nonlinear coupling constant 
u. While in Sections (II) and (III) we are dealing with the 
evolution of an initial state with fixed momentum, cho- 
sen at n = 0, i.e. ■!/'„(0) — (5(0), Section (IV) is devoted 
to the effects of a finite spread of initial conditions as 
strongly suggested by state-of-the-art experiments using 
ultracold atoms |ll[l9l|. 

The quantity we typically compute is the width of 
the momentum distribution of the system (p^(t)) — 



Y^'n=-oo '^^IV'nP, which gives the spreading of the wave 
packet over unperturbed levels or, equivalently, - apart 
from a constant factor 2 - the expectation value of 
the energy. The time-averaged spreading P^(r) = 

^ X]t=i {P^i^)) of the second moment has been frequently 
used (see e.g. 0|), as it preserves the exponent of the 
power law growth, while smoothing out oscillations. 

As pointed out in ^13.] the nonlinear shift is essential in 
determining dynamical features (providing for instance 
a mechanism for delocalization of the generic, irrational 
case) , but when we deal with delocalized states it is typ- 
ically quite small {e.g. for (M2) the shift is proportional 
to u/ An, where An is the width of the distribution over 
unperturbed states). So in general we may expect that, 
for moderate nonlinearities, the precise form of the shift 
does not alter in an essential way the nature of asymp- 
totic motion. 



II. LOCALIZED REGIME 

In this Section, we consider the regime where the value 
of the period of the (5-kick perturbation is incommensu- 
rate with 47r. For better comparison we fix t — 1 for all 
numerical computations shown in the following. The cor- 
respondent system in the u limit is characterized by 
the phenomenon of dynamical localization 16] caused by 
quantum interference effects. Previous theoretical pre- 
dictions and numerical simulations |]ji|20| indicate that, 
above a critical border Uc ~ 27r for the nonlinear cou- 
pling constant, dynamical localization is destroyed. The 
delocalization takes place in the form of anomalous sub- 
diffusion with an exponent of 2/5: in an asymptotic 
law (p^(t)) ^ c(u) • t'^^^ , (where c{u) ~ u'^^^) is predicted 
for both models; this is confirmed for both models by the 
data reported in Fig.^ In Fig.Qfa) a bilogarithmic plot 
of the time-averaged second moment vs time is shown for 
increasing values of the nonlinear coupling constant u. In 
spite of large oscillations, both modcls fit the predicted 
asymptotic behavior with a power-law exponent equal 
to 2/5. For nonlinear coupling larger than the critical 
border Uc, marked by the vertical line, the dependence 
of c(u) vs u is confirmed for both models by Fig. Q^b) . 
log{c{u)) is obtained by a one-parameter linear fitting of 
the logarithm of the second moment, once an anomalous 
diffusion exponent equal to 2/5 is assumed. 

The results obtained by calculating the evolution of the 
system with either (Ml) or (M2), starting from the same 
state ■0n(O) = (5(0) and parameters, appear to be different 
on short time scales; nevertheless the two models share 
the same asymptotic behavior of the time evolution of 
the second moment and of the dependence of c(u) vs u. 
The effect of the nonlinearity is the same for both mod- 
els only ed t — 1, because of the common initial state. 
As explained in the introduction, the way nonlinearity 
acts on the wave function is essentially different for (Ml) 
and (M2), at least before the state becomes delocalized: 
so deviations are qualitatively expected for intermediate 
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FIG. 1: (Color online) (a) Bilogarthmic plot of the time av- 
eraged second moment vs time in the localized regime. Time 
is measured in number of periods. The dashed and full lines 
refer to the (M2), with = 2^^, and (Ml), with = 2^^ and 
L = 80000, respectively. Values of u = 8, 10, 12, 14, 16, 20 
are considered; generally higher nonlinearity values yield big- 
ger spreading. The dashed line has the theoretically predicted 
slope 2/5. The values of the parameters are r = 1 and k = 2.5; 
the initial state is V'n(O) ~ S{0). (b) The logarithm of the co- 
efficient of the sub-diffusion as a function of log(u), for the sec- 
ond moment (circles) and its time-average (squares). Empty 
and full symbols refer respectively to (Ml) and (M2). The 
dashed lines show the predicted dependence ~ u"^^^. 

times, while we expect a closer analogy in the models' 
behavior in the asymptotic regime. Actually the close 
behavior exhibited by both models in Fig. Rafter a few 
time steps extends to more general features than the sec- 
ond moment: in Fig. |21we provide a comparison between 
full distributions over momentum states for u = 10. 

We remark that nonlinearity-induced delocalization 
has recently been explored also in stud ying survival prob- 
ability on a finite momentum sample |2ll |: while the au- 
thors use (M2) to cope with computational difficulties, 
our findings suggest that their results are probably rel- 
evant for true nonlinear Schrodinger equation dynamics 
too. 



III. RESONANT REGIME 

In this Section, we examine in detail the response of 
the system to the nonlinear perturbation in the resonant 
regime of the KR (t = Anr/q with r, q relatively prime 
integers) , characterized by a parabolic growth in time of 
the variance of the momentum distribution IT^ . The 
value of r is chosen equal to An, corresponding to the 



FIG. 2; Comparison between the momentum distributions for 
(Ml), circles, and (M2), full line, after t ^ W {&), t = 100 
(b), t = 1000 (c) kicks. The parameters are the same as in 
Fig-0 the nonlinear coupling is fixed to the value it — 10. 



first fundamental quantum resonance of the KR. 

In Fig. Ota) a bilogarithmic plot of the time-averaged 
second moment of the momentum distribution for 
different values of the nonlinearity is shown. The 
nonlinear coupling constant u varies from 1 to 400. As 
already noticed ,20i | , the resonant behavior survives even 
in the presence of nonlinearity, although generically 
the spreading is slowed w.r.t. the linear case. On 
asymptotically long time scales, the resonant growth 
with quadratic exponent is reached even for strong 
nonlinear perturbations, though we observe that the 
time needed to reach the asymptotic regime grows with 
u. The results in Fig. |3fa) obtained from (M2), shown 
in full lines, allow us now a more detailed analysis of 
the behavior of the system on quite long times. It can 
also be seen in Fig. ^a.) that the time evolution of the 
second moment of (Ml), shown with circle symbols, 
approaches the same asymptotic growth, even if some 
differences between the two models appear expecially for 
large nonlinearity (u > 50). 

The persistence of the resonant behavior in the pres- 
ence of nonlinearity can be explained intuitively as fol- 
lows. In the linear (i.e. u = 0) resonant case the width 
of the momentum distribution increases linearly in time. 
Therefore, from the normalization condition, the proba- 
bility amplitude to find the system in a momentum eigen- 
state n, decays as |V)„|2 ^ l/An ^ l/{nkt) 0,[l3. The 
nonlinear phase shift Tu\'ipn\'^ decreases with the same 
rate and its effects become irrelevant on long time scales, 
i.e. t » ^r. 
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FIG. 3: (Color online) (a) Resonant growth of the time aver- 
aged second moment vs time in the presence of nonlinearity, 
for T = An and k — 2.5. The initial momentum distribution 
is ?^„(0) = (5(0). The symbols and the full lines refer to (M2), 
with iV = 2", and (Ml), with N = 2^° and L = 5000000, re- 
spectively. The straight black line shows the resonant asymp- 
totic behavior . The values of the nonlinear parameter are 
M = 1, 5, 10, 20, 50, 100, 400. The inset is a magnifica- 
tion for u = 1 and m = 50 (lower part). A slight deviation 
between the two models can be seen for u — 50. (b) The 
function A'"'(t) vs time for (M2). Starting from above in the 
low t < 5000 region, the values of the nonlinear parameter 
are u = 0.1, 5, 10, 20, 30, 40, 50, 60, 70, 80, 100, 400. 



Nevertheless, the nonlinearity affects the evolution 
of the second moment on smaller time scales t < 
and introduces a w-dependence in the prefac- 
tor of the paraboHc growth law, {p^{t)) ^ a{u) ■ . 
Increasing nonlinear couplings manifests in a slower 
quadratic growth. In Fig. |nib) the function A*^")(i) — 
{{p^it)Y^^ - (p2(t))(°>) /t^ is plotted for different values 
of the nonlinearity. Increasing asymptotic absolute val- 
ues of A(tmax) give an estimate of the modifications in 
the transport induced by the nonlinearity. The coeffi- 
cients a A, calculated from the function A(t = tmax) with 
tmax = 20000 are shown in the inset of Fig. 0] 

A detailed analysis of the dependence of the coefficient 
a{u) of (p^it)) is presented in Fig. 01 The numerical 
calculation of log(a(u)) is obtained by a one-parameter 
linear fitting of the logarithm of the second moment vs 
the logarithm of time with a straight line of fixed slope 2. 
The fitting is performed on a time interval At = 200; this 
rather small time interval was chosen in order to make a 
comparison between the results from both models (empty 
and full symbols refer to (Ml) and (M2), respectively). 
The accordance between the two models is satisfying up 
to u < 50. For M > 50 the lowest order split method 
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FIG. 4: (Color online) Bilogarithmic plot of the coefficient 
a{u) of the quadratic growth of the second moment (circles) 
and its time-average (squares) as a function of the nonlinear 
parameter u. Empty and full symbols refer to (Ml) and (M2), 
respectively. The dashed lines show the algebraic behavior of 
a{u) for large u with an exponent equal to —4/5. Notice 
that for u ^ 1 the coefficients approach the theoretical values 
of the KR-model, marked by arrows. In the inset the empty 
symbols refer to a a calculated using the function A(t = tmax) 
with t^ax = 20000. 



[Tl| to evaluate the Floquet operator in (Ml) becomes 
less stable and the numerical errors around the borders 
of the finite basis propagate faster. 

Numerical data are compatible with an algebraic law 
a{u) ~ /c^/(2(l -I- u/c)''), where 7 is 4/5 and c is a con- 
stant of the order of 10 (for the time-averaged moment 
the constant fc^/2 is substituted by fc^/6). This law has 
the required asymptotic behavior for m — > 0: in this limit 
a{u) tends to the well-known value of the coefficient of 
the resonant KR, i.e. a(0) = P/2 . In Fig. HI the val- 
ues of a(0) are marked by arrows. For large values of u, 
a(u) decreases for increasing nonlinearity with the inverse 
power law u~^^^. At the moment we have no explana- 
tion for the minimum observed in the intermediate region 
{log{u) - 0.85). 

Up to now, we discussed only the case of attractive in- 
teractions, i.e. M > 0. It turns out that the fundamental 
quantum resonance at r = 47r is insensitive to the sign of 
the nonlinearity as it can be seen in the inset of Fig. [S] 
The same is true for the momentum distributions which 
are not presented here. On the other hand, the next or- 
der resonance at t = Stt is sensitive to the sign of u. 
For M < in Fig. I^b), the momentum distribution is 
slightly different from u > 0. Asymptotically, however, 
the same ballistic growth of the mean square momentum 
is obtained. This means that the details of the effect of 
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nonlinearity depend on the resonance type as far as the 
sign of u is concerned. This originates from the fact that 
while at the fundamental quantum resonances r = 47rrn 
(m > integer) the free evolution phase in the linear 
rotor is exactly one, at higher order resonances there is a 
nontrivial phase evolution between two successive kicks. 
The inset of Fig. IHIa) highlights that (p^)/2 either de- 
creases or increases with respect to the case u = at the 
second and fourth kick. This is related to the fact that 
at T = Stt the rephasing in momentum space occurs only 
every second kick, not between two successsive kicks as 
at T = Att. 

In Figs. (|5I6() the time scales relevant for experiments 
{t < 500) are investigated (results refer to (M2)). In Fig. 
(jSJ, we generically observe three regimes: (i) there is an 
initial stage, where the mean square momentum increases 
much slower than in the case u — 0. This stage is fol- 
lowed by stage (ii) where the increase can be faster than 
ballistic, and the mean square momentum can even be 
larger for larger nonlinearity (cf. \u\ = 10 as compared 
to \u\ = 1 in Fig.EJ. 

The observed superballistic growth of the second moment 
of the momentum distribution is quite surprising, in par- 
ticular, having in mind that such a growth is forbidden 
in the usual KR (i.e. u = 0) The results in Fig. 

are reminiscent of the observed superballistic spreading 
in ID tight-binding models however, here the su- 
perballistic behavior is caused by the nonlinear term in 
the time evolution, in contrast to the linear Hamiltonian 
models in . In terms of the model studied in 23] , the 
nonlinearity u would act as finite size trapping region 
(cf. also |2j|), outside of which the motion is ballistic 
(we already showed how nonlinearity does not essentially 
modify high n components). 

The final stage (iii) we call asymptotic regime, because 
there the growth exponent approaches the one for van- 
ishing nonlinearity (only for u = —100 this stage is not 
yet reached in Fig. EJ. 

In Fig. EJb) we notice the distinct peaks close to the 
very edge of the momentum distribution, for u = 0.2. 
Such peaks have been found for sufficiently large kicking 
strength fc > 2.5 and corresponding |u(fc)|=0.2...2, and 
it turns out that they can be up to one order of magnitude 
higher that the maximum of the momentum distribution 
for the linear KR at the resonances t = An and r — Stt. 
Fig. [71 compares the momentum distribution IV'nl^ at the 
fundamental quantum resonance r = An, and at the res- 
onance T = Stt for small nonlinearity u — —0.2 with the 
case of the linear KR. The distributions are shown after 
50 and 200 kicks, respectively, to stress their evolution 
in time. For both resonances, we observe a very interest- 
ing feature, namely, the small nonlinearity sharpens the 
edge peaks, which move ballistically, i.e. with a speed 
which is proportional to the number of kicks t (we recall 
that when u = the distribution is characterized by a 
largest momentum component also moving according to 
a linear law Umaxit) ~ ktn/2 Q). The peaks are more 
pronounced than in the linear case, and are remarkably 
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FIG. 5: Kinetic energy at the fundamental quantum reso- 
nance r = An, k — 7r,pi„itiai = 0, and nonlinearities u — 
(solid), —0.2 (diamonds), —1 (dashed), —10 (dotted), and 
— 100 (dash-dotted). The short dashed line shows the su- 
perballistic increase oc t'^-^ for the case u = —10. The inset 
presents the results for u — 0.2, 1, 10, and the case of m = —10 
(thin solid line) for better comparison. 
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FIG. 6: (a) kinetic energy as a function of the number of kicks, 
and (b) the corresponding momentum distributions after 200 
kicks, for the quantum resonance t — 3n, same parameters as 
in Fig. 1^1 apart the nonlinearity which is -S = —0.2 (diamonds) 
and 0.2 ((a) dotted, (b) solid). In (a) we show also data for 
u — (dashed) for comparison, and the inset illustrates the 
opposite effect of the second and fourth kick (the rephasing in 
momentum space occurs for the u = case only every second 
kick at T = Stt) depending on the sign of u. 



stable, i.e. their height decreases very slowly with in- 
creasing number of kicks in Fig. EJc), or even increases 
initially as in Fig. Oa). While we focussed our discus- 
sion on the model (M2) , the structure of the probability 
distribution is quite similar for (Ml), see Fig. (|SJ). 

The intermediate time scaling properties look in this 
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FIG. 7: Momentum distributions for zero initial momentum, 
k = TV, u — —0.2 (a,c) and u = (b,d), and t = Att (a,b), 
T = 3n (c,d). Note the stable peaks at the largest momenta in 
the case with small nonlinearity. The distribution are shown 
after 50 (dotted) and 200 kicks (solid line) in each panel. 





FIG. 8: (a) momentum distribution after t = 50 kicks, and 
(b) mean square momentum for pinitiai — 0, u — —0.2 and 
r = 4n, k = 77, shown are the evolutions induced by (M2) as 
diamonds, and by (Ml) as full lines, respectively. 



case more complex than in the kicked rotator dynamics 
in the presence of sticking accelerator modes, where the 
same exponent appears both in the classical and in the 
quantum case (where a new modulation appears). Work 
is in progress to see whether there exist classical map- 
pings that reproduce the peak dynamics we observe in 
intermediate time quantum behavior p3 |. 



IV. MOMENTUM DIFFUSION IN THE 
RESONANT REGIME 

All the above results have been obtained for an initial 
state in the form of a plane wave of null momentum. In 
a typical experiment, one can create a Bose condensate 
with an initial spread of momentum which is much less 
than two photon recoils (which can be imparted as mo- 
mentum kicks to the atoms by the kicking laser). We 
ask ourselves what happens if such a spread is taken into 
account. The momentum variable p, varing on a discrete 
lattice in the case of single rotor, becomes a continuous 
variable. For the linear KR, the eigenvalues of p can 
be written, distinguishing the integer and fractional part 
(quasi- momentum), as p — [p] + {p} = n + (3. Owing 
to the conservation of quasi-momentum (3, the system 
dynamics can be decomposed in a bundle of rotors 
(called in the following /3-rotors), each parametrized by a 
value of the quasi-momentum, evolving incoherently with 
operators with the same functional form of jSJl and 
in which h is substituted by n + p. Such a decomposition 
is not easily accomplished when we introduce a nonlinear 
term in the dynamics. The general task we have to face 
becomes the study of a nonlinear evolution equation with 
periodic coefficients. This is a quite a complex problem 
that cannot be tackled in full generality, even though dif- 
ferent approximation schemes have been proposed, e.g. 
by mapping the problem into a discrete lattice, which 
turns out to be useful if the wave function is expanded in 
a suitable set of localized functions related to the linear 
problem We generalize (M2) in such a way that its 
linear limit is the evolution operator corresponding to a 
quasimomentum /3 (as formerly specified), and assume 
that each nonlinear (3 rotor evolves independently. In 
this way we study the influence of nonlinearity on realis- 
tic initial conditions in a highly simplified way, by means 
of a generalized (M2) model: further work is obviously 
needed to check whether our findings extend to a full GP 
dynamics. 

The quantum resonance phenomenon in the KR is 
strongly sensitive to the values of the parameters of the 
system. The linear KR rotor (it = 0) exhibits the quan- 
tum resonance only for a finite set of quasi-momenta, i.e. 
[3 = ^ m/2p with m < 2p 0,[ll|2i|. A slight devia- 
tion of the quasi-momentum from f}^ changes completely 
the evolution of the system. For values of /? 7^ /3^, after a 
transient regime, the suppression of the resonant growth 
of the energy of the linear /3-rotor through the dynamical 
localization occurs; at fixed time t, only quasi-momenta 
within an interval ~ \/t oi (3^ mimic the ballistic behav- 
ior oc t^ and a rough estimate of the time up to which 
the quadratic growth of the /3-rotor energy persists is 
i ~ l/A/3, where A/3 = |/3 - /3^| [lEI. 

In the following, we investigate the mean square mo- 
mentum distribution of the generalized (M2) model, first 
for a) for a single (3 rotor with fixed quasi-momentum, 
and then b) for an incoherent ensemble of /3-rotors whose 
initial state in momentum space is a Gaussian distribu- 
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FIG. 9; (Color online) (a) Bilogarithmic plot of the second 
moment of a single /3-rotor with a fixed quasi-momentum for 
increasing nonlinear coupling (from below, referring to high 
t values, u — 0,1, 5, 10, 50, 100). The two lines have slope 2 
and 2/5. Note that in the KR case {u — 0), the localization 
occurs for time t> 1//3. The parameters are r — Att, k — 2.5 
and A'^ = 2^^. The initial state is an eigenfunction of the 
momentum with n — and /3 = 0.00946556. In the inset the 
calculations are prolonged ten times. The unbounded growth 
for u — 50, 100 can be clearly seen, (b) The same as (a) with 
u fixed {u = 5) and variable (3 (from above, referring to high 
t values, /3 = 0.0001,0.001,0.01 and 0.1). The arrows mark 
the times 1//3. 



(M2). The behavior in the correspondent hnear case of 
the KR is theoretically known (see Appendix A of [2^): 
for u = 0, the kinetic energy of the system increases diffu- 
sively in time with a coefficient proportional to and 
dependent on the initial distribution of quasi-momenta 
0. The presence of the nonlinearity manifests itself in 
a faster than linear growth, at least on short time inter- 
vals. After this transient regime, the asymptotic growth 
is expected to become approximately linear. The black 
straight line is drawn for better comparison. At fixed 
time t and assuming an uniform distribution of the quasi- 
momenta, the resonant rotors, whose quasi-momenta lie 
within the interval A/3, enter in the average of (p^)/? with 
a contribution of ~ t" and a weigth w ~ 1/t, while the 
non-resonant rotors give a contribution of (l~w)*t'^ . The 
exponents of the transport in the limit t —>■ +oo reach the 
values a{oo) = 2 and 2/5^7^1. Therefore, asymp- 
totically in time, the global transport exponent reaches 
the value 1. In the inset of Fig. ^| the exponents of the 
algebraic growth of the second moment are plotted as a 
function of the nonlinear coupling constant it. The fitting 
time intervals is 1000 kicks. Full and open circles refer to 
5000 and 500 quasi-momenta of the initial Gaussian dis- 
tribution: a slight rise in the exponents can be noted in- 
creasing the number of quasi-momenta, because a greater 
amount between them approches the value (3^ — 0, yield- 
ing the quadratic growth of the /3-rotor energy. Note that 
the exponent approches faster the value 1 for a uniform 
initial distribution of quasi-momenta (stars), confirming 
the previous argument. 



tion with zero mean and rms spreading A/3 — a — 0.01. 
We choose t = An, and the resonance condition is then 
met for /3^ = and /3^ = 1/2. The case with /3^ = 
fixed was considered in Section III. As in the localized 
regime, considered in Section II, the introduction of the 
nonlinearity causes a delocalization in the system with 
a non-resonant value of the quasi-momentum (/3 P^)- 
For small nonlinearities, the appearence of an anomalous 
asymptotic diffusion with an exponent of 2/5, after the 
initial ballistic behavior, is confirmed by data of Fig.Efa) 
for a /3-rotor with /3 sa 0.009. On the contrary, greater 
nonlinear couplings {u 50) introduce an excitation of 
diffusive type, starting from the first kicks. In Fig.|5{b) 
the quasi-momentum of the /3-rotor is varied and u is 
kept fixed. The arrows mark the times t, depending on 
the value of /3, approximately bounding the region of the 
ballistic growth. 

We then consider the dynamics of an incoherent en- 
semble of /3-rotors. The mean square displacement of 
the distribution is (p^(i))/3 = JdP(pl{t)). The average 
over /3 has been calculated using 5000 quasi-momenta. 
In Fig. ^1 the time evolution of the averaged second mo- 
ment of the initially Gaussian wave packet is shown for 



Fig. ^Jpresents a closer look at dynamics on a shorter 
time scale: the results refer to the case in which we found 
stable momentum peaks in Fig. [7| Part of the peak is 
still preserved for the used spread A/3 ~ 0.01, which can 
be realized in state-of-the-art experiments |li, ilSj . After 
about 15 kicks, more weight lies, however, now in the cen- 
tre of the distribution made up of rotors which do not ex- 
actly fulfill rephasing condition due to nonzero quasimo- 
menta. Also the increase of the mean square momentum, 
which is averaged incoherently over all the independently 
evolved initial conditions, is then not any more quadratic 
but closer to linear (see inset in Fig. Ill|l . as was found in 
the case of a uniform initial distribution of quasimomenta 
for the u = case 0, 113 • The mean square momentum 
still increases much faster than for nonresonant values 
of the kicking period r, where dynamical localization oc- 
curs. The latter may be destroyed by the nonlinearity but 
the above observed growth of (p^)/2 oc i^/^ (cf. Section 
II) is much slower than linear. On short time scales thus 
quantum resonance is very robust with respect to nonlin- 
ear perturbations. If our incoherent superposition model 
is correct after some initial stage, the ballistic motion 
should cease but the dynamics will show the influence of 
the ballistic quantum resonant transport. 
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FIG. 10: (Color online) Average over 5000 quasi-momenta 
of the second moment of the distribution vs time for (M2); 
r = 47r and k — 2.5. Different values of the nonlinear param- 
eter are shown in different colours; starting from below (refer- 
ring to high t values) u = 200, 100, 60, 40, 20, 0, 1, 5, 10; 
the initial wave packet in Fourier space is a Gaussian distri- 
bution centred in n = with rms a = 0.01. The inset shows 
the power-law exponents of the second moment as a function 
of u. The fitting is performed on time intervals At = 1000 
(circles) and At — 6000 (squares). Open symbols refer to 
500 quasi-momenta. Stars refer to a uniform distribution of 
quasimomenta. 
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FIG. 11: Momentum distribution after t = 10 (solid), 30 
(dashed), and 50 kicks (dash-dotted), for the same parame- 
ters as in Fig. |7^a), but incoherently averaged over indepen- 
dently evolved initial conditions (Gaussian initial momentum 
distribution with rms a = 0.01 centred around n = 0). The 
inset shows the corresponding (p^)/3/2 as a function of the 
number of kicks for u = —0.2 (solid) and u — (dashed). 



V. CONCLUSIONS 



In summary, we have numerically analysed in great 
detail the quantum transport occuring in two nonlinear 
generalizations of the famous (5-kicked rotor model, with 
a cubic nonlinearity as present in the Gross-Pitaevskii 
equation. We confirm previous results in the regime of lo- 
calized transport, and show the validity of the predictions 
of for a wide range of nonlinear coupling strengths. 
In addition, we found that the quantum resonances of 
the kicked rotor are very stable with respect to the non- 
linear phase perturbation, which looses its effect in the 
asymptotic limit of large interaction times with the pe- 
riodic driving. Surprising phenomena like pronounced 
peaks in the momentum distributions at quantum res- 
onance, and superballistic intermittend growth of the 
mean square momentum have been found. Both phe- 
nomena are caused by the cubic nonlinearity in the evo- 
lution, which shows that the analysed models bear a rich 
dynamical behavior in parameter space. 



The experimental works on the kicked rotor using a 
Bose-Einstein condensate ij mostly concentrated on the 
short time behavior at quantum resonance or on the so- 
called anti-resonance, where the motion is exactly peri- 
odic in the case u = 0. But an experimental observa- 
tion of the ballistic quantum resonance dynamics up to 
10 . . .30 kicks seems possible, for small enough kicking 
strength k such as to avoid a too fast spread in momen- 
tum space which cannot be monitored by standard time- 
of- flight detection 0, |1| . Our results are fully consistent 
with the few published experimental data, which show 
that both resonance and antiresonant dynamics essen- 
tially survive the presents of small nonlinearities, apart 
from other effects which damp, for instance, the periodic 
oscillations at the anti-resonance. Such effects are e.g. 
the uncertainty of the centre of the initial momentum 
distribution, and fluctuations in the experimental kick- 
ing strength . 
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